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_ Abstract: By use of the potential reconstruction approach we obtain a series of asymp- 

, ' totically AdS (aAdS) black hole solutions in an Einstein-Maxwell-Dilaton (EMD) system. 

Basing on the solutions of the system, we reconstruct a semi-analytical holographic QCD 
(hQCD) model with a quadratic term in warped factor. We discuss some aspects of the 
hQCD model, in particular we calculate the free energy of two static color sources (a 
heavy quark-antiquark pair) which is an important order parameter to describe confine- 
ment /deconfinement phase transition. The behavior of the free energy with respect to 
temperature and chemical potential is studied. We find that in the hQCD model the 
deconfinement phase transition can be realized and a critical point occurs. The resulting 
phase diagram in the temperature-chemical potential T— [i plane is in quite good agreement 
with the one from recent lattice results and effective models of QCD. 
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1 Introduction 

Understanding strongly coupling QCD related to the RHIC and LHC experiments is at- 
tracting much attention. Although a powerful method for this subject, the lattice QCD, is 
being developed, when it comes to the famous sign problem, lattice techniques are not well 
adapted to the case with finite real chemical potential \i and no many results have been 
produced so far. The other effective approaches, models and effective field theories also 
suffer from some problems in the dense matter case. On the other hand, the gauge/gravity 
duality [1-4] developed in string theory, can offer new insights to hadron theory [5-12] 
and strongly interacting quark gluon plasma (QGP) [13-20] in top-down and bottom-up 
setups. The duality approach can easily deal with the dense matter problem at least for 
the case in the deconfinement phase. However, for the hadron phase, the current status 
is still not very satisfied because the builded phase diagram so far [2 1-24] is different from 
that of the real QCD [25-27]. One is expecting to build up a holographic model which 
can completely recover the phase diagram of the real QCD. 

If one considers a system with graviton and dilaton [28-32], one can obtain solutions 
of Einstein equations by using the potential reconstruction method proposed in [33, 34]. 
However, in previous works, the back reacted geometry of the hadronic phase is not fully 
identified: if one wants to consider the degrees of freedom of quarks, one should take U(l) 
gauge field into consideration in the bulk [21, 22, 35-41]. Adding the additional field to 
the graviton-dilaton system, the U(l) gauge field is dual to the baryon number current 
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Jd = ip^ {x)ip(x), one may generate a chemical potential by turning on an appropriate 
electric field in the black hole geometry. On the other hand, in QCD the effect of finite 
quark density is introduced by adding the term Jjj = \i ip^(x)ip(x) to the Lagrangian in the 
generating functional, so that the chemical potential [i appears as the source of the quark 
density operator. According to AdS /CFT correspondence, the source of a QCD operator 
in the generating functional is the boundary value of a dual field in the bulk; therefore, 
the chemical potential can be considered as the boundary value of the time component of 
a U(l) gauge field Am dual to the vector quark current. Note that the £7(1) denotes the 
gauge symmetry in the bulk and global symmetry on the boundary. The symmetry is held 
by physical requirement of preservation of baryon number. Motivated from the holographic 
description of unquenched QCD system, we will design a graviton-dilaton-U(l) gauge field 
system to accommodate the degrees of freedom in full QCD. 

The confinement /deconfinement phase transition in QCD phase diagram is a very in- 
teresting issue to play with holography. In holographic QCD, it has been widely believed 
that the confinement phase in the pure Yang-Mills theory corresponds to the AdS DA soli- 
ton in gravity and the deconfinement phase corresponds to the black D4 brane pointed out 
by [42] . [43] argued that deconfinement in hard wall and soft wall models occurs via a first 
order Hawking-Page type phase transition between a low temperature thermal AdS space 
and a high temperature black hole. [44] extended this discussion by studying (charged) 
AdS black holes with spherical or negative constant curvature horizon. The authors of [42] 
and [45] investigated the deconfinement phase transition by introducing hard wall in the 
AdS/Reissner-Nordstrom black-hole with using Hawking-Page phase transition. Recently, 
the authors of [46] and [47] carefully considered the correspondence between phases and 
gravity backgrounds proposed by [42] and argued that the alternative gravitational config- 
uration named localized soliton would be properly related to the deconfinement phase. The 
deconfinement transition can be realized as a Gregory-Laflamme type transition. There are 
also various studies on Polyakov loop [48-52] to investigate QCD phase structure within 
hQCD models. We will propose a hQCD model and study its phase structure by studying 
the behavior of Polyakov loop operator in this model. 

In this paper, we will extend the potential reconstruction approach [33, 34] to an 
Einstein-Maxwell-Dilaton (EMD) system and build up a hQCD model with a positive 
quadratic term in warped factor of the bulk metric. Furthermore, we will calculate the 
free energy of two connecting Polyakov loop operators from holographic point of view 
to check whether the holographic model can realize the deconfinement phase transition. 
We introduce the quadratic correction term in warped factor to construct the gravity 
configuration with asymptotic AdS UV behavior. The motivation of choosing the quadratic 
correction is that the quadratic correction plays important roles to realize various facts in 

2 2 

low energy QCD. The work [53] by introducing the dilaton with form of e~ c z can realize 
the Regge behavior of hadron spectrum which can not be achieved in hard wall hQCD 
model[5][54]. Andreev and Zakharov introduced a positive quadratic correction, e cz with 
z the holographic coordinate and c > 0, to the warp factor of AdSs geometry. It turns 
out it is helpful to realize the linear behavior in heavy quark potential [55]. The linear 
heavy quark potential can also be obtained by introducing other deformed warp factors, 
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e.g. the deformed warp factor which mimics the QCD running coupling [56], and the 
logarithmic correction with an explicit IR cutoff log z '^~ z [57] . To produce the linear 
Regge behavior of the hadron excitations, Karch, Katz, Son and Stephanov [53] proposed 
the soft-wall AdSs model or KKSS model by introducing a quadratic correction to dilaton 
background in the 5D meson action, whose effect in some sense looks like introducing a 

2 

negative quadratic correction, e~ cz in the warp factor of the AdSs geometry. However, 
it is worth mentioning that the model with a quadratic correction in the warp factor 
of the metric is not equivalent to the model with a quadratic correction in the dilaton 
background, although both models give the same effective potential of hadron spectrum in 
the IR region [28, 29] . A positive quadratic correction e cz in the dilaton background of the 
5D hadron action has also been used to investigate hadron spectra [58, 59], however, higher 
spin excitations in this background will lead to imaginary mass [12, 60] and they also are 
inconsistent in the gluon sector. There are the flavor models [61-63] taking into account 
the chiral condensation. In [64, 65], the modified five-dimensional metric at the infrared 
region is constructed to obtain a nontrivial dilaton solution, which incorporates the chiral 
symmetry breaking and linear confinement. The modified factor is leading terms 1 + k 2 z 2 

i 2 2 

in e . They do predict the mass spectra of resonance states in the pseudoscalar, scalar, 
vector and axial-vector mesons, which agree with experiment data. [66] also confirmed 
that the quadratic term in dilaton plays an important role in linear confinement in the 
meson sector. Recently, basd on the AdS Reissner-Nordstrom black hole metric the work 
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[67] introduced the e z term in the dilaton, it turns out to be helpful to realize the heavy 
quark potential. All these studies mentioned above do not consider effects of back reaction 
of dilation and/or a modified warped factor . The work [57] studied the topic in a non- 

1,2 2 

critical string framework with back reaction effects and found that the quadratic term e z 
in warped factor is helpful to achieve the Cornell potential. The work [33] further confirmed 
the conclusion from thermal hQCD perspective. In some sense, a positive quadratic term 
captures some important QCD features, although we still do not well understand. Therefore 
we choose the warped factor with a positive quadratic correction to build up an effective 
hQCD model. 

The organization of the paper is as follows. In Section 2 we briefly introduce the poten- 
tial reconstruction approach and apply it to an Einstein-Maxwell-Dilaton system (EMD). 
In Section 3, we impose asymptotical AdS boundary conditions and regularity requirements 
on generic black hole solutions and figure out general formulas of thermodynamic quantities 
of the black holes. In Sections 4 , we make use of the reconstruction approach to propose our 
hQCD model with a positive quadratic k 2 z 2 term in warped factor and calculate relevant 
thermodynamic quantities such as temperature, entropy, etc, of the background solution. 
We also simply check the stability of dilaton potential from AdS/CFT perspective. In 
Section 5, we investigate free energy of a heavy quark pair and vacuum expetation value 
(VEV) of Polyakov loop which is order parameter to describe the deconfinement phase 
transition. In this hQCD model, the deconfinement phase transition can be realized and a 
critical point appears. Section 6 is devoted to conclusions and discussions. 
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2 Einstein-Maxwell-Dilaton system 

Let us begin with the following 5D Einstein-Maxwell-Dilaton (EMD) action in string frame 
J d 5 x^s e - 2 * (r s + 4d^<f> - V S (<P) - Jj-e^F^F^ , (2.1) 



where G5 and g g are 5D Newtonian constant and effective gauge coupling constant, g and 
Vs(4>) are the 5D metric and dilaton potential in the string frame, respectively. 

To solve this system, it is convenient to transform this action to Einstein frame. If the 
metric in Einstein frame g^ u and its corresponding one in string frame g^ v are connected 
by the scaling transformation 
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9% = effg^, (2.2) 



one can derive the exact relation between two actions in string frame and Einstein frame 
j V^e- 2 * (r s + 4d^</> - V S {4>) - -^efF^F^ 



up to a total derivative term, where 



4„ ,„„ , „ , 1 



R E - - V E (<j>) - -^F^ 



(2.3) 



-4<A 

V s = V E e—. (2.4) 



Thus we have the action in Einstein frame 

S 5D = — ^ J d 5 x^z (r - - VeW - ^2 V") > ( 2 - 5 ) 

where F^ u = d^A u — d u A^ is Maxwell field. Note that here we do not consider the coupling 
between gauge field and dilaton field in Einstein frame. 

For our aim, we are looking for black hole solutions of the system. In string frame we 
suppose the metric is of the form 

ds\ = (-f(z)dt 2 + J^ + dx'dx'j , (2.6) 

with L the radius of AdSs. We will set G5 = 1 and g g L = 1 in Section 3, 4 and 5. In 
Einstein frame, the metric becomes 

da% = — - 2 (-f(z)dt 2 + J^j + dx'dx 1 ) . (2.7) 

In Einstein frame, the gravitational field equations read 

E^u + \g% (U^t + V E (<t>)) - \d^d v 4> - ^ (f^fI - ^F P(T fA = 0, 

(2.8) 
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where E^ v is Einstein tensor. In the metric (2.7), the (t, t), (z,z) and (x\,xi) components 
of the gravitational field equations are respectively 

b'(z)f(z) 2 , 2 b{zfV E ^{z)) A' t (zf 

b {z) + ~wr + 9 b(z)(p {z) + em + = °' (2 - 9) 

2 _ M{zf _ _ 36(z) 2 VeM*)) 3A' t (zf 

9 1 J 4b(z)f(z) 4f(z) 16glb(z)2f(z) ' 1 ' ' 

6b'(z)f'(z) 6f(z)b"(z) 4 2 2 ^(z) 2 n 

/ (*) + 6(z) + + 3/ W (*) + &(*) WW) ~ = 0, 

(2.11) 

where b(z) = L2e z2 E , Ae(^) = j4 s (z) — \<j>(z) and .At (2) is electrical potential of Maxwell 
field. As a consistent check, turning off At{z) in Eqs. (2.9) -(2.11), one can easily reproduce 
equations of motion about the graviton-dilaton system discussed in [33]. 

Note that the above three equations are not independent. In (2.9)-(2.11), there are 
only two independent functions, therefore one of the three equations can be used to check 
the consistence of solutions. From those three equations one can obtain following two 
equations which do not concern the dilaton potential Ve(4>), 

ah, s t ,, ,f^'(z) 2\ 26"{z) U'(z) 

A" s (z) + A' s (z) I -*±± + - J - A' s (zf - -^J. - = (2.12) 



f"(z) + f{z) [SA'M - 24>'(z) -ty Z -± ggL2 = °- ( 2 - 13 ) 

Eq.(2.12) is our starting point to find exact solutions of the system. Note that Eq.(2.12) 
in the EMD system is the same as the one in the Einstein-dilaton system considered in 
[33, 34] and the last term in Eq.(2.13) is an additional contribution from electrical field. 
The equation of motion of the dilaton field is given as 

( jZ 3A s (z)-2<f> f ( \ \ L 5 5A a (z)-f 

ld z ^ 6 z3 !{Z) d z ^ - 9<pVE = o. (2 . 14 ) 

And the equation of motion for the Maxwell field is 

d.V^F^ = 0. (2.15) 



2 



E 



It is similar the equation given by [41]. Note that in (2.9)-(2.11) we have only considered 
the gravitational configurations with electric charge. 

We can see from the equations of motion of the system that once given a geometric 
structure A s (z), one can derive a generic solution to the system in Einstein frame. The 



- 5 - 



generic solution takes the form 



z e 2Mx) (lfo X y 2 e- 2A ^A' s (yfdy + ( i )1 ) 3A(z> 
y ax H \- 4>q, (2.16) 



o 



A t (z) = A t0 + A tl I J ye 2J ^- A ^ dyj , (2.17) 

/(*) = I xV^H^M 1 + fl \ dx + /„, (2.18) 

H>(.») 2A (z) 

V E (z) = 6 3 L2 3 (r 2 /" (*) - 4/(z) (3z 2 A';(z) - 2z 2 <fi"(z) + zV (z) 2 + 3) 



3zie i£M_ 2Mz)A , {z) , 



(2.19) 

where 4>q, (fii, A t o, A t \, fo, fi are some integration constants, and will be determined by 
suitable UV and IR boundary conditions. Generally, one cannot give the explicit form of 
Ve ((/>)■ But for some special cases, we can use the generating function A s (z) to obtain some 
analytical solutions of the graviton-dilaton-electric field system. In appendix, we give two 
analytical solutions which are generated by this potential reconstruction approach. 

For a consistent check, one can reproduce the general solution of graviton-dilaton 
system given by [33] by turning off the electrical field At(z). An alternative way to check is 
that one can set A s = 0, (fi = 0, f{z) = l,A t (z) = to obtain a constant dilaton potential 
Ve(z) = — as expected. 

3 Black hole solutions and associated thermodynamics 

In this section, we will work out general formulas of some thermodynamical quantities for 
the semi-analytical gravity solutions of the graviton-dilaton-electrical field system by using 
Eq.(2.16)-(2.19) for given metric ansatz in Einstein frame (2.7). Here we are interested in 
a series of solutions whose UV behavior is asymptotic AdS§ (aAdS). We also impose the 
requirements: /(0) = 1, and cfi(z), f(z), A t (z) are regular from z = to z^. Here Zh is 
supposed to be the black hole horizon with f(zh) = 0. In addition, we impose A t (zh) = 0, 
which is due to the physical requirement: A^A^ = g tt AtAt must be finite at the black hole 
horizon z = Zh- 

The first is how to parameterize the Hawking temperature of the black hole solution, 
which is defined by J £ ' . A black hole solution with a regular horizon is characterized by 
the existence of a surface z = Zh, where f(zh) = 0. The Euclidean version of the solution 
is defined only for < z < z^, in order to avoid the conical singularity, the periodicity of 
the Euclidean time is required as 

This determines the temperature of the solution as 

4-7T 
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With help of /(0) = 1 and black hole horizon z^, there f(zh) = 0, we can fix the 
integration constants Jq and /i in metric function f(z) in eq.(2.18) as follows. 

A 2 tl Jo 9(x) [Jo 11 g(r)dr g(y)Uy^j dx J* g(x)dx 
f{z) ~ l + ^p J* g(x)dx fi h 9(x)dx' (3 ' 3) 

where /(0) = 1, /i = - ft 1 ' ^ te g (x)fi g (y)is d y +i ; and the f unct i n <?(x) is defined as 

9fl9 Jo g(x)ax 

g(x) = x 3 e 20(x)-3A s (x)^ (34) 

One can easily check that f(zh) = from Eq.(3.3). One should confirm that there is no 
other Zh satisfying f(zh) = in the region < z < Zh- From Eq. (3.2), one can easily read 
out the relation between temperature and the black hole horizon from Eq.(3.3) as 



A n 9(zh) fo h g{r)dr f* h gl (y)dy - g(z h ) 



4vr 5 2L2 f Zh g(x)dx 



(3.5) 



Following the standard Bekenstein-Hawking entropy formula [68], from the geometry 
given in Eq.(2.7), we obtain the black hole entropy density s, which is defined by the area 
Aarea of the horizon 

A r, 3 



4G 5 V3 4G 5 




(3.6) 



where V3 is the volume of the black hole spatial directions spanned by coordinates Xj in 
(2.7). Note that the entropy density is determined in terms of horizon area in Einstein 
frame. 

In order to find out exact expressions about chemical potential and charge density, we 
should impose proper boundary conditions on At. Expanding At near z = 0, we have 

A t (z) = A t0 + A*e*P"W> (l + y^ - A'(y))) z* + ... (3.7) 

V 3 / y=0 

Having considered the boundary condition that At(zh) = 0, we can obtain from Eq.(3.7) 
the chemical potential \x and the integration constant Ati, which is related to the black 
hole charge, 



A t0 = \x (3i 
A tl = ^ = f^—. (3.9) 



Clearly the integration constant An can be determined by chemical potential \i and horizon 
Zh- From the coefficient of the z 2 term in (3.7), one can obtain the charge density of the 
black hole configuration. 

Before ending this section, we would like to stress that the integration constants ap- 
pearing in cf)(z) and Ve(z) do not occur in the solutions of electrical field A t (z) and f(z). 
As stated in [34], some integration constants will make contribution to the dilaton potential 
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Ve (</>), but this does not affect application of the potential reconstruction approach to the 
EMD system. Based on the builded gravitational configurations, one can build up various 
hQCD models by choosing proper A s . In this paper, we only focus on construction of 
holographic model to realize the deconfinement phase transition by investigating Polyakov 
loop. In addition, to be of physical interest, the builded balck hole solutions should be 
thermodynamical and dynamical stable. Therefore it is important to check the stability of 
the black hole solutions generally. In the next section, we will show that the mass of scalar 
field in our hQCD model indeed satisfies the Breitenlohner-Freedman (BF) bound. 

4 The hQCD model with a quadratic correction in warped factor 

From various works of constructing holographic QCD models for describing the heavy quark 
potential and the light hadron spectra, we learn that a quadratic background correction 
is related to the confinement property, i.e. the linear quark anti-quark potential [57] and 
the linear Regge behavior [12]. A positive quadratic correction, e cz with c > 0, in the 
deformed warp factor of AdSs can help to realize the linear heavy quark potential [55]. 
A quadratic dilaton background in the 5D meson action, whose effect in some sense looks 

2 

like introducing a negative quadratic correction, e~ cz , in the warp factor of the AdSs 
geometry, is helpful to realize the linear Regge behavior of hadron excitations [53]. In our 
previous study in [33], it shows that the hQCD model with a positive quadratic correction 
in warped factor is favored. Therefore, we introduce the following hQCD model with a 
positive quadratic correction in the warp factor of AdSs i n Eq.(2.6), i.e. we take 



where k is a parameter related to energy scale of the model, which will be fixed later. A s {z) 
is the generating function of the solution (2.16)-(2.19). In terms of the general expressions 
for the solution of graviton-dilaton-electric field system, we can obtain the configuration of 
4> field as 



where we have set (fro = to keep the metrics (2.6) and (2.7) are asymptotical AdS, and 
H(z) is 



A s (z) = k 2 z 2 



(4.1) 



d>(z) = + H{z)) 



(4.2) 




(4.3) 



With A s and </>, we can obtain the metric function / as 



f( z ) =f{z,z h ,n) 



9l L2 \f* h g(y)Uy 
Jo g(x)dx 
Io h 9{x)dx ' 



) 



2 Jo 9( x ) [lo h 9(r)dr f* g(y) I dy^j dx 
Io h 9(x)dx 



(4.4) 
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where 



5 (x)=x 3 e§* 2 * 2 ( 1 +"(*))- 3fe2 * 2 . (4.5) 
With the solutions of A s {z) and (j)(z), one can obtain the expression of the electric field 

A t ( z ) = fj,+ ^— f xe**^- 1+H Wdx. (4.6) 

Further one can obtain the potential of the dilaton field. The potential includes the contri- 
bution from the electrical field. The form of the potential is too complicated and therefore 
we do not present it here. Instead we will discuss the UV behavior of the dilaton potential 
to check whether the potential satisfies the constraint from the 5D Breitenlohner-Freedman 
(BF) bound. 

The conformal invariance in the UV can be restored when (f> ~ at the UV boundary 
z — > 0. One can expand (/>(z) at UV boundary z ~ as 

i / ~\ 3 A* z 3k z , A 

<t>{* -»> 0) — + + ■■■ ■ (4-7) 

The behavior shown in Eq.(4.7) is consistent with the requirement of the asymptotic AdSs 
near the ultraviolet boundary. 

Through the AdS/CFT dictionary, for any scalar field we have 

Urn V($) = -§ + ^ A(A - 4)* 2 + 0($ 4 ). (4.8) 
By using the following relationship 

«V(.)-**r*»^> , )-Mj + 0(.) + ..., (4.9) 



d<& dr \ <9<3> dr 

one can easily get the conformal dimension of the scalar field. In Eq. (4.8), A is defined 
as A(A — 4) = M^L 2 , which is constrained by the BF bound 2 < A < 4. Comparing our 

case to the standard formula, one should notice the realtion = \ %6, thus we have 



Ml = - A (4.10) 

Therefore, the conformal dimension of the dilaton field is A = 2 in our model. The dilaton 
therefore satisfies the BF bound but does not correspond to any local, gauge invariant 
operator in 4D QCD. Although there have been some discussions in recent years of the 
possible relevance of a dimension two condensate in the form of a gluon mass term [69], 
it is not clear whether we can associate 4> with dimension-2 gluon condensate, because the 
AdS/CFT correspondence requires that bulk fields should be dual to gauge- invariant local 
operators. This situation is the same as the case in the recent paper [33]. Note that in the 
potential reconstruction approach, the electrical field does not change the configuration of 
dilaton and the UV behavior of the dilaton potential. 
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Once the metric function / is given, one can easily find the relation between the black 
hole temperature and the horizon radius Zh- In figure 1 we plot the relation between the 
temperature and horizon with different chemical potential. One can see from figure 1 that 
the black hole temperature can not reach zero temperature if the chemical potential is 
zero. On the other hand, if one turns on a nonvanishing chemical potential, the lowest 
temperature can reach zero, this case corresponds to an extremal black hole. In the case 
with a nonvanishing chemical potential, as one can see from the figure that the behavior of 
the temperature heavily depends on the value of chemical potential. Basically the behavior 
of the temperature can be classified into two cases: when chemical potential fj, > fj, c , the 
temperature decreases monotonically to zero as z^ increases, while if fi < fj, c , it decreases 
monotonically to a minimum at z m , then goes to a maximum at zm and then decreases 
to zero. The critical chemical potential depends on the parameter k. In our model, we 
take parameter k = 0.3GeV, which is consistent with the lattice data for equation of state, 
heavy quark potential and so on [33]. We have checked that if one changes the value of 
k, our conclusions in this paper do not change qualitatively. In the case of k = 0.3GeV, 
the corresponding critical chemical potential fj, c = 0.34GeV. Here the critical chemical 
potential is determined by T(/i = fj, c , Zh = z m ) = T(/i = /i c , zh = zm)- And z m and zm 
have been presented in figure 2. 



T(GeV) 




Figure 1. The black hole temperature versus horizon Zh with different chemical potential. When 
fx = 0, the temperature behavior is similar to the case considered in[33]. When [i > fi c the 
temperature monotonically decreases to zero as Zh increases. While < li < fx c , the temperature 
decreases to a minimum at z m and then grows up to a maximum at zm and then decreases to zero 
monotonically. When fi = fi c , one has z m — zm- z m and zm are marked explicitly in figure 2. In 
our model, we fix k = 0.3GeV and accordingly the critical chemical potential is [i c = 0.340GeV. 

From the behavior of the temperature, one can see that the black hole solution is 
locally thermodynamical stable when fj, > fj, c , because in this case, the heat capacity of the 
solution is positive (Note that the horizon size of the black hole is \ j Zh). On the other hand, 
the case with \jl < \i c has to be discussed separately. To clearly see the stability of the black 
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hole solution in the case \i < fx c , in figure 2 we plot the temperature versus the horizon z^ 
in the case fi = O.lGeV < fi c . We can see from the figure that the background black hole 
is thermodynamical unstable in the region z m < < zm, where z m and zm are the black 
hole horizons, respectively, corresponding to the minimal and maximal temperatures. In 
this region, the heat capacity of the black hole is negative. The black hole solutions in the 
regions Zh < z m and Zf- L > zm are thermodynamical stable. When ji = fj, c , z m and zm are 
degenerated to one point. 

T/GeV 




Figure 2. The temperature of the black hole with /i = O.lGeV. The three black hole solutions 
with horizon z' p , z p and z' p ' have the same temperature. The black hole with z m < z p < zm is 
thcrmodynamically unstable. In the figure we take g g L = 1, k = 0.3GcV. 

The existence of the critical chemical potential plays a crucial role to realize the critical 
point in QCD phase diagram. We will see this in the next section by calculating the heavy 
quark potential in this hQCD model. 

5 Heavy quark potential, Polyakov loop and QCD Phase diagram 

To investigate some properties of the hQCD model constructed in the previous section, 
let us study an infinitely heavy quark-antiquark pair, at distance r from each other. It is 
interesting to investigate how the free energy of such a system changes with temperature 
and chemical potential by using the holographic description of loop operators [70-72]. 
The free energy is a proper quantity to describe the confinement /deconfinement phase 
transition. 

5.1 Heavy quark potential and Polyakov loop 

At finite temperature, the free energy F(r, T) of an infinitely heavy quark-antiquark pair 
at distance r can be obtained in QCD from the correlation function of two Polyakov loops 
[72] 

(P(f 1 )pt(x 2 )) =e -^^T) +7 (T) (51) 
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with r = \x\ — X2\ and 7(T) a normalization constant. Moreover, the vacuum expectation 
value of a single Polyakov loop 

(V) = e -^ F °°^ (5.2) 

(F°°(T) = F(r = oo,T) and the normalization factor is neglected) is an order parameter 
for the deconfinement transition of a gauge theory [72]. Within the gauge/string duality 
approach, we can attempt a calculation of the expectation values in (5.1) and (5.2) by 
considering a fundamental string configuration with Polyakov loop on the boundary. The 
Nambu-Goto action of the string is 

Sng= 2^J d ^^tH = / d 2 ^det[g MN (d a X^) (d b X*)] (5.3) 

In terms of the geometry background gMN given by (2.6), the induced metric g a b can be 
expressed as 



9ab = p 2A s (*)(ld£ , A , (5.4) 




where we have chosen the static gauge £° = t and £ l = x. We are interested in the 
configurations of two static quarks on the boundary. The configurations should obey the 
boundary conditions: z(x = 0) = zq and z'{x = 0) = and z(x = ±|) = 0. Here the prime 
stands for the derivative with respect to x. There are two independent configurations as 
shown in figure 3 which satisfy the boundary conditions. Figure 3(a) stands for the case 
with two static quarks linked by a fundamental string. This configuration describes the 
confining phase. Figure 3(b) describes the case where the string is broken into two segments, 
each of them extends to the horizon of black hole. This configuration stands for two static 
decoupled quarks, which corresponds to the deconfinement phase. The configuration in 
figure (a) corresponds to the "Minkowski embedding", while figure (b) to the "black hole 
embedding" in the D3/D7 setup [73]. 

From (5.3) and (5.4), we obtain the free energy of the string configuration 



F(r,T) = 2* f dx^- y/f(z) + (zf , (5.5) 
vr J_ r/2 z z v 



f0 n k 2 z 
<■/ 

L- 



where g p = — r. One can see from (5.5) that the configuration of the string satisfies the 
constraint 

* = ,5,) 
^Jf(z) + (z'f 

where H is a conservation quantity, the Hamiltonian of the string configuration. 

For later convenience, we express F(r,T) in terms of zq and /o = f(z)\ z=zo . Defining 
v = j^, and subtracting the UV (v — > 0) divergence in (5.5), which corresponds to sub- 
tract the infinite quark and antiquark mass in four dimensional QCD [70], we obtain the 
regularized free energy for the string configuration in the confining phase 
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Figure 3. The two configurations represent the confinement phase and deconfincmcnt phase. In 
Figure (a) there is a fundmental string connecting with two static quarks on the boundary. This 
configuration denotes the confining phase. In Figure (b) one string beaks into two segments, both of 
them extend to the horizon of the black hole. This configuration corresponds to the deconfinement 
phase. 



IT A 

where F = F/k, A = k zq and 

r{v) = J 



. 1 dv 
-1+ / — 

o v 



a X 2 v 2 



t(v) 



fo 



/ f(zov) 

The distance r = kr can also be parameterized by zq, 



4 e 2X 2 (l-v 2 ) 



r(z ) 



2X 2 Vfo 



dv 



2 p X 2 (l-v 2 ) 



v e 



(5.7) 



(5.8) 



(5.9) 



'o r(v)f(z v) 

Next we calculate the regularized free energy for the configuration in the decconfine- 
ment phase. In this case, one has zq = Zh- Using the gauge £° = t and £ x = z and x'(z) = 0, 
we obtain the regularized free energy 



F°°(f) = — 

7T 



1 



2ft 



dz 



— + / — (e 

Zh 



+ C(A,r), 



(5.10) 



where T = T/k and z = kz. T) is a term related to the regularization procedure, and 
we identify the maximum of the free energy F(r, T) with F°° to fix the regularization term 
following [67]. 

5.2 QCD phase diagram 

In this subsection, we will discuss the phase diagram of the hQCD model by analyzing 
the regularized free energies obtained in the previous subsection for those two different 
configurations. 
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To analyze the behavior of f(zo) is helpful to distinguish two different phases. Note 
that there might exist a divergence in (5.9) for some parameters due to a vanishing t(v). 
An infinite f in figure 2(a) implies that two static quark are still coupled even with infinite 
separation. This means that the case with the infinite f is in the confining phase. If there 
does not exist any divergence in r(zo), one can find a maximum r max at some zq with 
fixed chemical potential \x and temperature T. Beyond f max , the string configuration in 
the confining phase becomes unstable, and the deconfinement phase transition happens in 
this case. 

We can see from (5.8) that when v = 1, one has r = 0. This means that the divergence 
in (5.9) appears when r = 0. Note that other factors in (5.9) do not lead to any divergence. 
Let us first analyze the behavior of t(v) near v = 1. Expanding the quantity in the square 
root in t(v) at v = 1, one has 



t(v) ~ Vci(l -v) + c 2 (l - v) 2 + 0((1 - v) 2 ), (5.11) 
where c\ and c 2 are two expansion coefficients. Note that the following relations 

1 , 1 

dv , = 2, 
o y/l-v 

1 1 

dv = oo. (5.12) 

o 1 -v 

We can see that if c\ = 0, the integration in (5.9)) is divergent, otherwise it always gives a 
finite result. Furthermore we find that c\ is of the following form 

z df(z ,z h ,fi) 2 2 

ci~— + 8kz -4, (5.13) 

J{zo,Zh,n) dz 

with positive Z/,,/z and < zq < z^. Note that the first term in the right hand side in 
(5.13) is negative due to d/O^-fe^) < in the region zq < Zh, and the second term is always 
positive. Therefore, in the case with fixed fx, if z^ is large enough, c\ can reach zero at 
some zo(< Zh)- At that point, t(zq) is divergent. On the other hand, in some region of 
parameters Zh and fj,, c\ is always positive definite and r(zo) is then finite in that region. 
Let us stress again that the case with a vanishing c\ is always in the confining phase, 
while the case with positive definite c\ corresponds to deconfinement phase. In figure 4 we 
show the behavior of c\ with respect to zq for various z^ and a fixed chemical potential 
fi = O.lGeV. 

One can also obtain free energy F of the static quark-antiquark potential in two differ- 
ent phases through (5.7). In figure 5, we plot the function r(zo) with a fixed [i = O.lGeV. 
We see clearly that when z^ < z p and then c± > 0, t(zq) has a maximum r max . This case 
corresponds to the deconfinement phase. On the other hand, the z^ > z p cases correspond 
to the confining configuration in which there exists a divergence in the blue curve v(zq). 
In the confining phase, if the separation t{zq) of two static quarks can go to infinity, the 
phase is called permanent confinement phase. The red and pink solid curves correspond 
to the deconfinement case. In these cases, there exist some maximums in the curves r(zo) 
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Figure 4. This plot shows the behavior of c\ with respect to zq with three different black hole 
horizons Zh and a fixed chemical potential /i = O.fGeV The solid red curve represents the decon- 
finement phase, the solid black curve corresponds to the confinement phase, and the solid blue 
one denotes the deconfinement phase transition. z p is the critical horizon which leads c\ to have 
a single zero root, zh < z p and zh > z p correspond to the deconfinement phase and confinement 
phase respectively. Here we take g g L = 1 and k = 0.3GcV. 




Figure 5. This plot shows the relation between the distance between two quarks and zq with a 
fixed (i = O.fGeV. The blue curve corresponds to the configuration given in figure 3 (a) which 
stands for confining phase, while the red and pink curves correspond to the configurations given 
figure 3 (b). When Zh = z p , t(zq) can be infinity at some zq in < zq < Zh- The dashed curves are 
related to configurations in unstable region. 

in < zq < Zh, beyond those maximums there are no stable configurations of two coupled 
quarks. 

In figure 6, we plot the free energy with respect to zq with a fixed fi = O.lGeV. We 
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can see from the figure that in the deconfinement case, the free energy has a maximum 
in < zq < Zfr. We identify the maximum of F(r,T) with F°°(T) following the strategy 
given in [67]. In figure 7 the free energy is plotted with respect to the distance between 
two quarks. We see from the figure that the potential with Zh = z p can go to infinity when 
r — > oo. This means the vacuum expectation value (VEV) of polyakov loop (5.2) in this 
phase is vanishing, which implies the configuration is in the confining phase. On the other 
hand, the free energy F in the cases Zh = 0.7z p and Zh = 0.9z p is finite when r — > oo. 
This implies that in this case the VEV of Polyakov loop (5.2) is not vanishing and the 
configuration is in the deconfinement phase. 

F(GeV) 

1.5 r 




Figure 6. The free energy F of string configurations versus zq with a fixed [i = 0.1 GeV. The solid 
blue curve corresponds to the confinement phase and the other two to the deconfinement phase. 
Here g g L = l,g p — l,k = 0.3GeV. The dashed curves are related to configurations in unstable 
region. 

Through the above analysis, we see that the key point of finding phase transition in 
the hQCD model is to find the critical horizon Zh = z p , for which t(zq) is divergent. The 
strategy is the same as the one employed in works [67, 74]. The numerical results in figure 
5 show that one can fix the critical horizon z p and then the phase transition temperature 
with a fixed chemical potential fj, through (5.13). 

Let us first discuss the case with fi < fi c . In this case the phase transition is a first order 
one. To confirm this, let us see the phase transition from the Polyakov loop perspective 
[75]. From our numerical study, one can see that the critical horizon z p determined by 
free energy of heavy quark pair always lives in z m < z p < Zm when \i < fi c . The phase 
transition temperature corresponds to the black hole having the horizon z p . But the black 
hole with z p is thermodynamical unstable. Furthermore one can see from figure 2 that 
there exist another two black hole solutions with horizon z' p and z p , which have the same 
temperature as the black hole with horizon z p . The two black hole solutions with zL and 
z'p are thermodynamical stable. They correspond to two different phases, deconfinement 
phase and confinement phase, respectively. The deconfinement phase transition happens 
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0.5 



0.0 
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-0.5 
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Figure 7. This plot show the free energy F(r, T) versus the distance r between two quarks with a 
fixed fi = 0.1. The solid blue curve corresponds to the confinement phase and the other two to the 
deconfinement phase. In the confining phase the free energy F goes to infinity when r — > oo. In 
the deconfining phase there is a maximum in < zq < Zh- Here g g L = l,g p = 1, k = 0.3GeV. The 
dashed curves are related to configurations in unstable region. 

from the black hole with z 1 ' to the one with z' p at the transition temperature. 

After fixing the phase transition temperature, next we calculate the VEV of a single 
Polyakov loop by using (5.10), which is a function of temperature. We find that when 
Zh > z p , (V) is always vanishing, while (V) does not as Zh < z p . Therefore there must 
be a jump of (V) from to a nonvanishing one in (V) — T plane as shown in figure 8(a). 
The black dot and red dot in the figure represent the black hole solutions with z' p ' and z' p , 
respectively. We can clearly see from the figure that the phase transition is a first order 
one in this case. 

Now we discuss the case with \i > \x c . In this case, the temperature is a monotonic 
function of horizon z^ as shown in figure 1. After fixing the phase transition temperature 
through (5.13), we calculate the VEV of the a single Polyakov loop by using (5.10). We find 
that when z^ > z p , (V) = 0, which corresponds in the confining phase, while {V) increases 
monotonically from as Zh < z p . The behavior of (V) is shown in 8(b). This case can be 
interpreted as a continuous transition [76]. The continuity of (V) at the phase transition 
point implies that the deconfinement transition might be a crossover in the heavy quark 
limit [77]. 

Combing the above analysis, we plot a naive T — fx phase diagram in figure 9 for 
our hQCD model. In the figure the black dot denotes the critical point with \i = \i c . 
When [i < fi c , the deconfinement transition is a first order one, while it is a continuous 
transition as /i > [i c . This phase diagram agrees with the expectation from effective 
field theory [78, 79] and recent lattice QCD simulation [76]. By using effective theory, 
the recent lattice QCD simulation [76] studies the deconfinement transition of QCD with 
heavy quark. It is found that the deconfinement transition is first order in small chemical 
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potential region, while it is an analytic crossover in large [i region. Our results therefore are 
consistent with the lattice simulation [76]. Furthermore, note that the work [67] considered 
a hQCD model by introducing a warped factor to deform the AdS Reinsser-Nordstrom black 
hole background and concluded that the deconfinement phase transition is a continuous 
transition. Differing from the work [67], we have constructed a self-consistent gravitational 
configuration by considering the back reaction of dilaton field and a critical point has been 
found in the phase diagram. 




(a) (b) 



Figure 8. The vacuum expectation of a single Polyakov loop versus temperature. The left panel 
shows the case with [i < [i c and the right one with /i > fi c . In (a) the dashed pink curve shows the 
behavior of (V) in the unstable region z m < Zh < zm • Here g p = 1, k = 0.3GeV. 




Figure 9. The phase diagram of the hQCD model constructed in this paper. The dashed blue 
curve stands for the continuous phase transition and the solid blue one for the first order phase 
transition. The black dot represents the critical point. Here g p = 1, k = 0.3GcV. 
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6 Conclusion and discussion 

The gravity /gauge theory duality is a powerful tool to investigate strongly coupling sys- 
tems. In this spirit, one of dreams is to build a holographic description dual to a real 
low energy QCD theory. In this paper, by using the potential reconstruction approach 
[33, 34], we have given a generic formulism to find a series of asymptotically AdS black 
hole solutions for the Einstein-Maxwell-dilaton system. In this approach, the back reaction 
of dialton and Maxwell field is taken into account. In this sense, our approach avoids some 
shortcomings in the literature in constructing holographic QCD models. 

Based on the approach, we have constructed a self-consistent gravitational configura- 
tion to describe some properties of low energy QCD theory. The gravitational configuration 
includes a quadratic term in warped factor of bulk metric. The behavior of temperature 
of the black hole configuration is similar to that of an AdS Reissner-Nordstrom black hole. 
The quadratic term in warped factor plays an important role in the hQCD model. By 
calculating heavy quark potential and Polyakov loop in this hQCD model, we have ana- 
lyzed the phase structure of the model. It has been found that there exists a critical point 
in T — fj> phase diagram. When the chemical potential \i < /i c , the deconfinement phase 
transition is a first order one, while it is a continuous transition when [i > fi c . This phase 
diagram is agreement with results from effective field theory [78, 79] and recent lattice 
QCD simulation [76]. In our model, the value of the critical chemical potential depends on 
the model parameter k. In our discussions, k = 0.3GeV and then [i c = 0.34GeV. We have 
checked other values of k and our results do not change qualitatively and the conclusions 
still holds. 

In this work we have only studied the deconfinement transition of the hQCD model in 
heavy quark limit, by studying Polyakov loop in the black hole background presented in 
this paper. It would be of great interest to further investigate other aspects of the hQCD 
model. For example, it is required in the model to further study the spectra of hadrons, 
chiral phase transition [38-40] , hydro dynamical properties of QGP, and color flavor locked 
phase [80-82], etc.. In addition, it would be also interesting to construct thermal gas 
solution in the EMD system and to discuss the Hawking-Page phase transition between 
the black hole solution and the thermal gas solution. 
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APPENDIX 



A Two analytical black hole solutions in EMD system 

In this appendix, we give two analytical black hole solutions of EMD system by using 
Eqs.(2.16)-(2.19). Here we are interested in the solutions whose UV behavior is asymptotic 
AdS§. We impose the constraint /(0) = 1 at z = 0, and require 4>(z) and f(z) to be 
regular in the region z £ [0, Zh], where z = is the AdS boundary and Zh corresponds to 
the horizon of black hole solution. 

We give black hole solutions in Einstein frame. The metric in Einstein frame takes the 
following form 

2 L 2 e 2As ~ a ( N , 9 dz" 2 



ds E = ~2 y-f(z)df + y^y + dx l dx l ) (A.l) 

T 2 2A E / 1 2 

-f(z)dt 2 + + dx l dx l ) , (A.2) 



z 2 v v 1 m 



where A E {z) = A s (z) - 



3 ' 

The first set of exact solutions is 



A E (z) = log 



z 



zosinh(^) 



4V U fn , ,z 4 z 6 1 2 . , ^ z 



4 

/(z) = i_^lli(3 s inh(— ) 4 + 2sinh(— ) e ) + -F^ sinh ' 

3 z zo 8 yzo 

A t (*) = /i - ^Vi 2 sinh , (A.3) 

where zq and [i are two integration constants, V\\ and V12 are two constants from the 
dilaton potential, and g g is gauge coupling. The dilaton potential is given by 

12 + 9 sinh 2 ( + lSVu sinh 6 ( t) V 2 sinh 6 ( ^■ 



Ve{4>) = ^ + 8L 2 • (A.4) 

If turn off the gauge field, one can reproduce one solution of the graviton-dilton system 
given in [33] [34]. On the other hand, if set V\\ = and V\2 = 0, one can reach the 5D BPS 
solution in [34]. 
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The second set of exact solutions is 

A E {z) = -l og (l + £) , (A.5) 
„ , , T/ / z 7 z 6 3z 5 z 4 \ 

p 2 zg / 5z 8 z 10 z 9 llz 7 -z 6 \ 
+ ^L 2 \32? Q + 604° + l^f + sZz 7 + 24zfJ' ( } 

0(z) = 3\/2sinh _ 

^(z)=/* + p(^p + y). (A.7) 

where zq, f-t, and p are integration constants and V21 is a constant from the dilaton potential 
and g g is gauge coupling. The dilaton potential is given as 



12 42sinh 4 (^) 42sinh 2 (^) 



L 2 L 2 L 2 

3V- 21 sinh^(^) 3y 21 sinh^(^) W 21 (^) 

35L 2 10L 2 10L 2 

5^ 



2 L 2 l 20L 2 + 10L 2 + 80L 2 



}• (Ai 



16L 2 8L 2 28L 2 

This solution is also a generalized one given in [33] . 

In this paper, we have not discussed thermodynamical properties of these two sets of 
black hole solutions and possible applications in hQCD model from gauge/gravity duality. 
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